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We analyze the impact of recent solar and atmospheric data in the determination of the neutrino oscillation 
parameters, taking into account that both the solar i/ e and the atmospheric may convert to a mixture of 
active and sterile neutrinos. Furthermore, in the context of the atmospheric neutrino problem we discuss an 
extended mechanism of neutrino propagation which combines both oscillations and non-standard neutrino-matter 
interactions. We use the most recent neutrino data, including the 1496-day Super-K solar and atmospheric data 
samples, the latest SNO spectral and day/night solar data, and the final MACRO atmospheric results. We confirm 
the clear preference of all the data for pure-active oscillation solutions, bounding the fraction of sterile neutrino 
involved in oscillations to be less than 52% in the solar sector and less than 40% in the atmospheric sector, 
at 3cr. For the atmospheric case we also derive a bound on the total amount of non-standard neutrino-matter 
interactions, bounding the flavor-changing component to —0.03 < e < 0.02 and the non-universal component to 
|e'| < 0.05. 



1. Introduction 

The experimental data on atmospheric neutri- 
nos jjJD show, in the muon-type events, a clear 
deficit which cannot be accounted for without in- 
voking non-standard neutrino physics. In addi- 
tion to this, the recent results from the Sudbury 
Neutrino Observatory (SNO) on neutral current 
(NC) events || have added more weight to the 
already robust evidence that an extension of the 
Standard Model of particle physics is necessary in 
the leptonic sector. Altogether, the most popular 
explanation of both the solar and the atmospheric 
neutrino problem is provided by the neutrino os- 
cillation hypothesis. However, many alternative 
attempts to account for neutrino anomalies with- 
out oscillations have recently been proposed j|. 
While the results of solar neutrino experiments 
still admit very good alternative explanations [^) , 
the atmospheric neutrino anomaly is so well re- 
produced by the — > v T oscillation hypothe- 
sis that one can use the robustness of this 
interpretation to place stringent limits on a num- 
ber of alternative mechanisms Q . 

In the present work we present an updated 
analysis of the solar and atmospheric neutrino 
data, first in the context of the standard oscil- 
lation hypothesis (Sees. @ and 0), and then as- 



suming that neutrino posses also non-standard 
interactions with matter (Sec. fl|). Motivated by 
the stringent limits from reactor experiments [^) , 
we adopt an effective two-neutrino approach in 
which solar and atmospheric analyses decouple. 
However, in the pure-oscillation case presented in 
Sees. U and || our effective two-neutrino oscillation 
approach is generalized in the sense that it takes 
into account that a light sterile neutrino, advo- 
cated to account for the LSND anomal y |9fl , may 
participate in the conversion process |1C|] . The 
natural setting for such a light sterile neutrino 
is provided by four- neutrino models. In this pa- 
per we will determine the constraints on the os- 
cillation parameters in this generalized scenario 
following from solar and atmospheric data sep- 
arately, addressing the reader to our recent pa- 
pers m Refs. for a mass-scheme-dependent 
combined analysis of all current oscillation data. 

2. Solar neutrino oscillations 

In the following we analyze solar neutrino data 
in the general framework of mixed active-sterile 
neutrino oscillations. In this case the electron 
neutrino produced in the sun converts into a com- 
bination of an active non-electron neutrino v a 
(which again is a combination of and v T ) and 
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Table 1 

Best fit values of Am 2 OL and ^ol, with the corresponding yJol an d the GOF for 81 — 2 d.o.f., for 
pure-active and pure-sterile solar neutrino oscillations. 



a sterile neutrino v s : 



\/\-T) s v a + y/r^Vs 



(1) 



The parameter rj s with < rj s < 1 describes the 
fraction of the sterile neutrino participating in the 
solar oscillations. Therefore, the oscillation prob- 
abilities depend on the three parameters Am 2 OL , 
9sol and r) s . 

As experimental data, we use the solar neu- 
trino rates of the chlorine experiment Homes- 
take jl3|], the most recent result of the gallium 
experiments SAGE @ and GALLEX/GNO (l5), 
as well as the 1496-days Super-Kamiokande data 
sample |16 in the form of 44 bins and the lat- 
est results from SNO presented in Ref . || , in the 
form of 34 data bins. Therefore, we have a total 
of 3 + 44 + 34 = 81 observables, which we fit in 
terms of the three parameters Am^, 6 soh and 
r) s . Further details about our statistical analysis 



can be found in |10|. 

We have performed a global fit of solar neutrino 
data, whose results are summarized in Tab. |l| 
and Fig. |l|. Our global best fit point occurs in 
the LMA region, for tan 2 SOL = 0.44, Am 2 OL = 
6.6 x 1CT 5 eV 2 and r) s = 0. We obtain a x r 2 nin = 
66.1, which for 81 — 3 d.o.f. leads to an excel- 
lent goodness of fit of 83%. From Fig. |l| we can 
derive the 3<r bounds 0.25 < tan 9 SOL < 0.83, 
2.6 x 10" 5 eV 2 < Ato 2 ol < 3.3 x 10" 4 eV 2 and 
r) s < 0.52. Note that maximal mixing c/sol — 45° 
is now ruled out, and that the upper bound on 
Am 2 OL is rather solid even without the inclusion 
of reactor experiments |8| . 

From Tab. we notice the strong discrimina- 
tion against non-LMA solutions implied by the 
present data: in the pure-active case we find 



Axlow = 9.0, A X 2 AC = 8.9, A X 2 MA = 23.2 and 
Ay 2 . „ 2 = 30.9 relative to the global best fit 

^Just-So ° 

in the LMA region. This shows that the first 
hints in favor of a globally preferred LMA oscil- 
lation solution, which followed mainly from the 
flatness of the Super-K spectra, have now become 
a robust result, thanks to the additional data to 
which SNO has contributed significantly. Note 
that especially SMA and Just-So 2 are highly dis- 
favored with respect to LMA. 

The inclusion of a sterile neutrino participating 
in the oscillation process is strongly disfavored by 
the solar data. From the right panel of Fig. [I] we 
first see how the preferred LMA solution survives 
in the presence of a small sterile component char- 
acterized by 77 s . However, increasing ij s leads to a 
clear deterioration of all the oscillation solutions. 
Notice that there is a crossing between the LMA 
and VAC solutions, as a result of which the best 
pure-sterile description lies in the vacuum regime. 
However, in the global analysis sterile oscillations 
with r) s = 1 are highly disfavored. We find a y 2 - 
difference between pure-active and pure-sterile of 
32.9 in the LMA region, which decreases to 22.9 
if we allow also for VAC. This means that pure- 
sterile oscillations are ruled out at 4.8a compared 
to the active case. 

In summary, we have found that, as long as 
the admixture of sterile neutrinos is acceptably 
small, the LMA is always the best of the oscilla- 
tion solutions, establishing its robustness also in 
our generalized oscillation scheme. 

3. Atmospheric neutrino oscillations 

In our analysis of atmospheric data we make 
use of the hierarchy Am 2 OL <C Am 2 TM and neglect 
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the solar mass splitting ||. Further, in order to 
comply with the strong constraints from reactor 
experiments Q we completely decouple the elec- 
tron neutrino from atmospheric oscillations Jl7[ . 
In the following we consider a generalized oscilla- 
tion scheme in which a light sterile neutrino takes 
part in the oscillations |l(J. Such a scenario re- 
quires two more parameters in addition to 9 ATU 
and Anij TU . We use the parameters and d s al- 
ready introduced in Ref. [[n] , and defined in such 
a way that (1 — d^) and (1 — d s ) correspond to the 
fractions of and v s participating in oscillations 
with Am 2 TM , respectively. Hence, pure-active at- 
mospheric oscillations are recovered in the limit 
dfj, = and d s = 1. 

In addition, we also present a "restricted" anal- 
ysis in which the is completely constrained to 
the atmospheric mass states, so that d^ = 0. In 
this limit the parameter d s has a similar interpre- 
tation as rj s introduced in the solar case, and 
oscillates into a linear combination of v r and v s : 

dfj, = : — > \fd s v T + \J 1 - d s v s . (2) 

For the atmospheric data analysis we use 
all the charged-current data from the Super 
Kamiokande |jj and MACRO Q experiments. 
The Super-Kamiokande data include the e-like 
and /^-like data samples of sub- and multi-GeV 
contained events (10 bins in zenith angle), as well 
as the stopping (5 angular bins) and through- 
going (10 angular bins) up-going muon data 



events. From MACRO we use the through- 
going muon sample divided in 10 angular bins 
Therefore, we have 65 observables, which we fit 
in terms of the four parameters A 

^ATM i ^ ATM i 

and d s . Further details can be found in Ref. [[LQL 
The results of our analysis are summarized in 
Fig. ||. In contrast to the solar case, the atmo- 
spheric x 2 exhibits a beautiful quadratic behav- 
ior, reflecting the high quality of the fit and the 
robustness of the oscillation solution. The global 
best fit point occurs at sin 2 9 ATU = 0.49, Am 2 TM = 
2.1 x 10~ 3 eV 2 , d s = 0.92 and d^ = 0.04, and ex- 
hibits a small but non-vanishing sterile neutrino 
component. However, this effect is not statisti- 
cally significant, since also the pure- active case 
with d s = 1 and = gives an excellent fit: 
for sin 2 ATM = 0.5 and Am 2 TM = 2.5 x 10~ 3 eV 2 , 
we have that the x 2 -difference with respect to the 
best fit point is only 3.3. 

From Fig. || we can extract the 3<r allowed 
ranges: 0.29 < sin 2 ATM < 0.71 and 9.2 x 
10~ 4 eV 2 < Am 2 TM < 4.6 x 10~ 3 eV 2 for the stan- 
dard two-neutrino oscillation parameters, while 
for d s and d^ we have (1 — d s ) < 0.40 and 
(1 — d tl ) > 0.85. So we see that atmospheric 
data essentially reject a sterile neutrino compo- 
nent, and strongly bound the ^-content in atmo- 
spheric oscillations to be nearly maximal. Note 
that if we impose the condition d^ = 0, the 3cr 
limit on the sterile neutrino fraction improves to 
(1 - d s ) < 0.19. 
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Figure 2. A%| TM as a function of Am ATM , sin 2 ATM , d s and d M . In each panel the undisplayed parameters 
are integrated out. The "Mixed" and "Restricted" cases correspond to d s = 0.5 and — 0, respectively. 



These limits on the sterile admixture are signifi- 
cantly stronger than obtained previously jflj , and 
play an important role in ruling out four-neutrino 
oscillation solutions in a combined global analysis 
of the LSND anomaly g§. 

4. Non-standard neutrino interactions 

In this section we still focus on atmospheric 
data, but instead of adding a sterile neutrino we 
now consider the possibility that neutrinos are 
massive and moreover possess non-standard inter- 
actions with matter (7) . This may be regarded as 
generic in a large class of theoretical models. For 
definiteness, in the following we assume that non- 
standard neutrino interactions occur only with 
the d-quark. Also, as in Sec. || we completely de- 
couple v e from atmospheric oscillations. In this 
case, the propagation of and v T inside the 
Earth is governed by the Hamiltonian: 

_ Ato 2 tm /- cos20 ATM sin20 ATM \ 

COS 2 $ ATM / 

±V2G F N d (r)(l £ £ ^j, 

where the sign + (— ) holds for neutrinos (anti- 
neutrinos), Nd(r) is the number density of the d- 
quark along the path r of the neutrinos propagat- 
ing in the Earth, and s and e' are phenomenologi- 
cal quantities describing flavor-changing and non- 
universal non-standard interactions, respectively. 



Therefore, as for the sterile neutrino case the 
transition mechanism depends on four indepen- 
dent parameters. Note that without loss of gen- 
erality we can restrict the range of the oscillation 
parameters to < ATM < 7r/4 and Am 2 TM > 0, 
provided that we consider both positive and neg- 
ative values of the NSI parameters e and e' 0. 
The details of our analysis are the same as in the 
previous section. 

Our results are summarized in Fig. ||. The 
global best fit point occurs at sin 2 (26* ATM ) = 1, 
Am 2 TM = 2.3 x 1(T 3 eV 2 , s = 6.7 x 1CT 3 and 
e' = ±l.lxl0~ 3 . As in the sterile neutrino case, 
a small component of NSI is preferred, however 
the effect is not statistically significant: the best 
pure oscillation solution e = e' = 0, which occurs 
at ATM = 45° and Am 2 TM = 2.5 x 10~ 3 eV 2 , ex- 
hibits a x 2 which is worse than the global one only 
by 2.4 units. Thus, the determination of the os- 
cillation parameters Am 2 TM and 9 ATM is stable un- 
der the perturbation introduced by the additional 
NSI mechanism, as can be seen also comparing 
the "Free" and "Oscillations" lines in the first 
two panels of Fig. [| Moreover, the \ 2 function 
is quite flat in the e' directions for e' —> 0, and 
almost symmetric under the exchange e' — > — e'. 

On the other hand, the pure NSI solution 
Am 2 TM = gives a very poor fit, and as already 
found in Ref. it is completely ruled out. This 
occurs since the NSI mechanism is not able to rec- 
oncile the anomaly observed in the upgoing muon 
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sample with that seen in the contained event sam- 
ple. Therefore, when combining the two mecha- 
nisms of — > v T transition, oscillations play the 
role of leading mechanism, while NSFs can only 
be present at a sub-dominant level. 

From Fig. || we also derive the 3cr allowed 
ranges: 0.84 < sin 2 (26» ATM ) < 1 and 1.0 x 
10~ 3 eV 2 < Am 2 TM < 4.8 x 10~ 3 eV 2 for the 
oscillation parameters, while for the NSI param- 
eters we have -0.03 < e < 0.02 and |e'| < 0.05. 
This is the main result of our analysis, since it 
provides limits to non-standard neutrino interac- 
tions which are truly model independent, being 
obtained from pure neutrino-physics processes. 
In particular they do not rely on any relation be- 
tween neutrinos and charged lepton interactions. 
Therefore our bounds are totally complementary 
to what may be derived on the basis of conven- 
tional accelerator experiments [jl8|| . 

5. Conclusions 

In this paper we have presented an updated 
analysis of solar and atmospheric data, both in 
the context of standard oscillation hypothesis 
(generalized to account also for a light sterile neu- 
trino), and assuming that neutrino posses also 
non-standard interactions with matter. 

We have found that both solar and atmospheric 
neutrino data are very well explained by the 
simplest hypothesis of oscillations into an active 



neutrino, and disfavor the inclusion of an extra 
sterile state. In fact, at the 3tr level the frac- 
tion of sterile neutrino involved in oscillations 
is limited to be less than 52% in the solar sec- 
tor and less than 40% in the atmospheric sec- 
tor. Furthermore, atmospheric data also con- 
strain non-standard neutrino interactions with 
matter, which are bounded to —0.03 < e < 0.02 
for the flavor-changing component and < 0.05 
for the non- universal component. 

In all the considered cases the determination 
of the oscillation parameters Am 2 and 9 is stable 
under the inclusion of exotic physics, such as an 
extra sterile neutrino or non-standard neutrino- 
matter interactions. This fact is an evidence 
in favor of the three-neutrino oscillation solution 
to the solar and atmospheric neutrino problems. 
In particular, for what concerns the atmospheric 
anomaly we conclude that a maximum mixing an- 
gle is a solid result, which must be incorporated 
into any acceptable particle physics model, even 
in the presence of exotic neutrino interactions. 
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